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Several anomalies have been identified which may imply a breakdown of the statistical isotropy of 
the cosmic microwave background (CMB). In particular, an anomalous alignment of the quadrupolc 
and octopole and a hemispherical power asymmetry have increased in significance as the data have 
improved. There have been several attempts to explain these observations which explore isotropy 
breaking mechanisms within the early universe, but little attention has been given to the possibility 
that these anomalies have their origin within the local universe. We explore such a mechanism by 
considering the kinetic Sunyaev-Zel'dovich effect due to a gaseous halo associated with the Milky 
Way. Considering several physical models of an anisotropic free electron optical depth contributed 
by such a halo, we find that the associated screening maps of the primordial anisotropics have the 
necessary orientations to affect the anomaly statistics very significantly, but only if the column 
density of free electrons in the halo is at least an order of magnitude higher than indicated by 
current observations. 



I. INTRODUCTION 

High-precision measurements of the cosmic microwave 
background (CMB) anisotropics and the agreement of 
these observations with our theoretical expectations has 
elevated cosmology to a precision science. The observed 
anisotropics have been measured to follow isotropic, 
Gaussian statistics with an almost scale-invariant power- 
spectrum. However, as observations have improved, 
anomalies which might indicate deviations from the the- 
oretical expectations have increased in significance. Pro- 
saic explanations for the observed anomalies may yet be 
found as our understanding of the systematics of the ob- 
servations improves, but it is also of interest to explore 
whether these anomalies point us towards exotic and un- 
expected physics. 

One such anomaly is the apparent breakdown of statis- 
tical isotropy that has been reported in the CMB fluctu- 
ations at the largest observable scales [IHS] measured by 
the Wilkinson Microwave Anisotropy Probe (WMAP). 
Two major observations 1 that are suggestive of this sta- 
tistical anisotropy refer to particular special directions 
in the sky. The first is the planarity of the quadrupolc 
and octopole and their mutual alignment; this plane is 
also roughly orthogonal to the CMB dipole direction 
[U [S] . The second is an asymmetry in the amplitude 
of the power spectrum between the northern and south- 
ern ecliptic hemispheres [3] which has also been observed 
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A quadrupolar anisotropy of the inferred primordial power spec- 
trum has also been reported 6-8;, which is aligned with the eclip- 
tic. However, this signal shows significant variations between 
different WMAP differencing assemblies at the same frequency 
[7], pointing to an experimental systematic which has not been 
taken into account in these analyses. 



in Cosmic Background Explorer (COBE) data but at a 
lower significance [3]. 

With the evidence for these anomalies increasing as the 
data improves [51 UOj many studies have proposed mod- 
ifications of early-universe physics in order to generate 
a violation of statistical isotropy (see Refs. [TTUI4] for a 
small selection) . There have also been discussions of how 
a violation of statistical isotropy would affect other obser- 
vations assuming that its source is primordial rT3HT7] , 
enabling consistency tests of those hypotheses. 

However, soon after these anomalies were discovered, 
it was noticed that the directions associated with them 
corresponded to special directions within the local uni- 
verse: four of the planes associated with the quadrupole 
and octopole are orthogonal to the ecliptic, the remaining 
octopole plane is orthogonal to the Galactic plane, and 
the hemispherical asymmetry is aligned with the ecliptic. 
The heuristic connection between these CMB anomalies 
and our local environment has generated only a handful 
of quantitative attempts to connect the two, including 
possible foregrounds associated with the heliosphere [18] , 
the kinetic Sunyaev-Zel'dovich effect of free electrons in 
the Galactic disk [19, 20), the thermal Sunyaev-Zel'dovich 
effect of the local universe 21, 22 , the Rees-Sciama effect 
of the local superclusters [23] . and the Integrated Sachs- 
Wolfe effect of the low-redshift universe [21] ■ This is an 
area which deserves greater attention, since a very local 
signal can affect the largest observable scales. 

In this paper, we explore another possible local ori- 
gin of the anomalies by investigating how the scattering 
of CMB photons by free electrons diffusely distributed 
within the Milky Way halo through the kinetic Sunyaev- 
Zel'dovich (kSZ) effect may introduce an anisotropic con- 
tamination of the primordial signal and help explain the 
origin of the anomalies and their alignment with special 
directions in the local universe. The effects of the kSZ 
from the Galactic halo on the CMB has been recently 
discussed in Ref. [25] . However Ref. [25] did not explore 
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how the kSZ may impact these anomalies and instead 
concentrated on whether the local kSZ may produce a 
foreground that must be subtracted in order to extract 
the primordial signal. 

The prediction of an extended gaseous halo within the 
Milky Way is quite robust. Models of galaxy formation 
for halos with masses M > 10 12 Af Q predict infalling 
gas is shock-heated to the virial temperature (around 
10 6 — 10 7 K for the Milky Way) and remains in hydro- 
static equilibrium until it is able to cool and condense 
to form stars [2"6l |2"7] . For a Milky Way-sized halo it is 
expected that the baryonic mass fraction follows the cos- 
mic average, / ~ 0.1. By subtracting the total baryonic 
mass in stars and gas in the Galaxy, a reasonable esti- 
mate of the baryonic mass in an extended, hot, gaseous 
halo is - 5 x 10 10 M [15]. Observations of OVII and 
OVIII X-ray absorption have also indicated the existence 
of an extended hot gaseous halo |28 associated with the 
Milky Way. Furthermore, observations of pulsars yield 
estimates for the column density of free electrons and 
lead to a free-electron fraction within the Milky Way halo 
which is consistent with a hot gaseous halo with a column 
density of ~ 10 21 cm" 2 [29]. Given that the Milky Way is 
moving relative to the CMB with a velocity v/c ~ 10~ 2 , 
the local kSZ may produce a signal around 1 /jK. 

In § [II] we point out that any anisotropic optical 
depth to Thomson scattering off local electrons, coupled 
through the kSZ effect with the dipole of the local elec- 
trons with respect to the CMB, induces an anisotropic 
imprint with a black-body spectrum. Its alignment is 
determined by the dipole direction and the anisotropy of 
the distribution. Given that the kSZ is expected to be 
subdominant to the cosmological signal, in § |III| we point 
out the disproportionate impact on anomaly statistics of 
small signals. We make no judgement on the usefulness 
or otherwise of such a posteriori anomaly statistics (for a 
post-WMAP 7 year discussion of this point, see Ref. 30 ). 
Instead, we take at face value anomaly statistics which 
have been invoked in a large body of literature, and study 
the impact of the kSZ signal on their significance. In § |IV| 
we derive the coupling of an anisotropic electron screen 
to the Doppler dipole, and consider several physical mod- 
els of the form of the anisotropic electron distribution in 
§ W\ We compute the impact of these halo kSZ models 
on anomaly statistics in § |VI| and discuss the results in 
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II. 



EFFECT OF LOCAL KSZ ON CMB 
ANISOTROPIES 



A simple way to understand how an anisotropic optical 
depth leads to additional anisotropy in the CMB is to 
note that along a given line of sight, n, with optical depth 
r(n), the fraction of photons which scatter out of the line 
of sight is given by 
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where AT(n) is the temperature anisotropy along that 
line of sight. In addition, there are photons that scatter 
in to the line of sight isotropically from every other part 
of the sky; thus, they contribute the mean temperature, 
T, and the observed temperature anisotropy is |31j 
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From this simple calculation we can see that an 
anisotropic optical depth will couple with any inherent 
temperature anisotropies and thereby modulate the ob- 
served anisotropy. 

A more rigorous way to derive the same result is to 
consider the Boltzmann equation which dictates the evo- 
lution of the photon distribution function. In cosmol- 
ogy the application of the Boltzmann equation leads to 
the usual evolution equations for the photon distribu- 
tion function within a perturbed Friedmann-Roberson- 
Walker (FRW) universe (32J, [33] . It can also be used to 
describe the evolution of the photon distribution func- 
tion as the photons pass through the Milky Way's halo. 
Denoting perturbations to the photon temperature by 
T®(t,x,n) the Boltzmann equation gives 
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where the first term describes scattering out of the line 
of sight, the second term is due to scattering into the 
line of sight, the last term is the Doppler shift due to 
the bulk motion of the electrons (with velocity Vb in the 
CMB rest- frame) . The Thomson scattering cross-section 
is ot = 6.65 x 10 - 25 cm 2 . 

Any optical depth contributed by free electrons within 
the Milky Way halo will be much smaller than unity so 
we may solve Eq. pi perturbatively. The zeroth order 
solution follows d0W /dt = 0. We then find the first 
order correction to this 
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where fj, = n • k and S^ ' is the term in parentheses 
in Eq. ([3| evaluated for the zeroth-order solution. In 
the case of the CMB this term is approximately given by 
5(0) = -n-v b since v b /c = 6x 10" 3 [M] and 6^ - 10~ 5 . 
As can be verified by substitution, the solution to Eq. Q 
is then given by 
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where the integral extends along the line of sight. The 
observed temperature pattern is then given by 
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where we have explicitly written 8'°' = 9 cos and n e is 
the number density of electrons within the Milky Way 
halo. Defining C(n) = <tt ds n e (sn), the amplitude 
of the second term in Eq. (|6|) can be written 
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where we have defined the angle averaged C = 
fd 2 nC{n)/(4n). 

Without any reference to a specific model for the op- 
tical depth, it is clear that the kSZ is capable of produc- 
ing any modulation of the CMB that we may observe. 
Letting r(n) denote the optical depth as a function of 
position on the sky assuming some underlying primor- 
dial temperature anisotropics then the optical depth is 
determined by 
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where D is the Wigner matrix corresponding to the ap- 
propriate rotation between the frames, and 4> can take 
any value. The statistics given by the maximum values 
of L\, L\ and = L\ + L\ have been used in the lit- 
erature to quantify the planarity of the quadrupole and 
octopole, and also to capture their mutual alignment, re- 
spectively. 

Table |T] shows the values and p- values for these statis- 
tics for the WMAP 7 year Internal Linear Combina- 
tion (ILC7) foreground-cleaned map |36 , with and with- 
out subtracting the kinematic quadrupole (KQ). The 
p-values are computed from 10,000 isotropic realiza- 
tions. The kinematic quadrupole induced by the dipole 
anisotropy due to our motion relative to the CMB is 
at the level of —1 /jK —> +3 fiK, whereas the CMB 
quadrupole is at the —22 fiK —> +16 /iK level. The Ta- 
ble shows that the subtraction of the very subdominant 
KQ signal leads to a factor of ~ 2.7 decrease in the p- 
value of the most anomalous statistic, L\ z . This effect is 
even more dramatic for the alternative WMAP 5 year 
Tegmark-Oliveira-Costa-Hamilton (TOH5) foreground- 
cleaned map [2J, where the maximum value of L\ z goes 
from 0.93 to 0.96 upon subtracting the KQ. The p-value 
corresponding to the latter is 0.0029, making the statistic 
a factor of 5.5 more unlikely under the null hypothesis of 
isotropy. This dramatic impact of subtracting the small 
KQ signal was also noticed in the WMAP first year data 
by Ref. [37], and is also present in the WMAP 5 year 
ILC (ILC5) map [38]. For purposes of comparison, the 
p-value for the ILC5-KQ L\ 3 = 0.95 is 0.0062. 



III. DISPROPORTIONATE IMPACT OF SMALL 
SIGNALS ON ANOMALY STATISTICS 

While, in principle, this effect can generate the ob- 
served large-angle isotropy anomalies, in practice (given 
the expected column density of free electrons and the 
dipole velocity of the CMB), estimates of the kSZ are 
at the 1 fxK level (compared to the ±20 ^tK of the 
quadrupole). Given that, at most, the kSZ can produce 
~ 10% modulations of the primary CMB, it would appear 
as though its effect on the observed anomalies would be 
negligible. However, as we now discuss, even a 10% mod- 
ulation can have a large effect on the inferred statistical 
significance of these anomalies. Consider the "angular 
momentum dispersion" statistic [5] [T31 |3S] for quantify- 
ing the planarity of multipole I: 
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which, given the observed realization on the sky, is maxi- 
mized in some direction ft . The transformation between 
the general frame and the maximizing frame is given by 
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TABLE I: Values and p-values (computed from 10,000 
isotropic realizations) for angular momentum dispersion 
statistics L|, L% and L| 3 (see text) for the WMAP 7 year 
Internal Linear Combination (ILC7) map, with or without 
subtracting the kinematic quadrupole (KQ). 



Statistic 


ILC7 


ILC7 


ILC7-KQ ILC7-KQ 




maximum 


p- value 


maximum p-value 


r 2 


0.943 


0.66 


0.983 0.38 




0.919 


0.17 


0.926 0.15 




0.931 


0.015 


0.953 0.0055 



Given these observations, it is very important to quan- 
tify the impact of the kSZ effect of the Galactic halo for 
two reasons. First, the velocity of the Galactic barycen- 
tcr with respect to the CMB points in the direction 
(l,b) = (266.5°, 29.1°) [39 , which is quite close to the 
Solar CMB dipole direction, (263.8°, 48.2°) and thus rel- 
evant for the orientations picked out by the large-angle 
isotropy anomalies (see Fig. [T]). Secondly, measurements 
of the column density of free electrons in the Galactic 
halo are not very precise and come primarily from theo- 
retical arguments; observational limits are very difficult 
to obtain and the uncertainties on constraints are large 
|40) . Thus, in the following calculations, we will consider 
a number of plausible physical models of the geometry of 
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the halo free electron optical depth, but allow the magni- 
tude of the effect to float somewhat, roughly at the level 
of the KQ amplitude, while keeping it greatly subdomi- 
nant to the CMB quadrupole/octopole amplitudes. 




FIG. 1: The WMAP 5 year quadrupole and octopole, with 
the ecliptic plane (solid line), the Solar system dipole (SD), 
and the Galactic barycenter dipole (GD). 



IV. COUPLING TO DOPPLER DIPOLE 

After processing through an anisotropic Thomson scat- 
tering screen, an anisotropic temperature field gets 
screened as 0(n) — > 0(n) exp(— r(n)), as we have seen. 
Consider that the electron distribution consists of an 
isotropic part f and a (small) anisotropy r scr (n). This 
can be expanded as e -1 "^ ~ e~ r (l — T scr (n)), lead- 
ing to an observed temperature anisotropy obs (n) = 
e~ r [1 — r SCI (n)] 0(n), with the screened component of 
the temperature field T scr (n) = e~ T r scr (n)0(n). 

Let us expand the anisotropic part of the local optical 
depth in spherical harmonics, 

r s "(h)= £$> imi y, imi (n), (11) 

£ t >l mi 

which we will couple to the Doppler dipole of the CMB, 
T di P(n) = A d Y w (n). (12) 



The screened temperature field after passing through the 
anisotropic r scr , assuming this is much smaller than the 
isotropic part of the screening optical depth, is then given 
by 
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Expanding T scr in spherical hamonics, we obtain 
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The right hand side of this equation can be expressed in 
terms of Wigner 3j symbols using the Gaunt integral 
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where the second expression has made use of the symme- 
try property of the Gaunt integral that m 1 +m 2 +m 3 = 0. 
The first 3j symbol in Eq. ( 16 ) enforces the triangle con- 



dition l\ = £ ± 1, so finally we obtain 
where 
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Remember that in these expressions, the Tg" m » starts at 

C = 1. So, Tom OC Cl m Tl m , T\ m (X Ci m Ti m , T2m oc 

(CtmTXm + C 3m r 3m ), and so on. 



V. MODELS FOR ANISOTROPIC FREE 
ELECTRON DISTRIBUTIONS 

In this section we consider how an extended gaseous 
partially ionized halo couples the CMB to the Doppler 



dipole through the kSZ effect, first gaining physical intu- 
ition using an analytic toy model of an anisotropic optical 
depth distribution, and then considering physical models. 



A. Simple geometries for the electron distribution 

Let us develop a feel for some simple, plausible elec- 
tron distributions. Consider t(6,4>), a smooth function 
defined on the unit sphere, with (0 < < it, < <fi < 2tt), 
which we want to expand in terms of spherical harmonics 
as 
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The spherical harmonic coefficients for r are given as 
usual by 
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A sphere only has the non-zero coefficient aoo so is not 
useful for our purpose. Consider the ellipsoid centered 
on the origin and aligned with the coordinate axes, 

where as usual, x = r sin 9 cos tf>, y = t sin 9 sin (f>, z = 
t cos 9. Thus for this distribution, 



= ac (a 2 cos 2 9 + c 2 sin 2 ( 
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This function only depends on 9, so the only non-zero 
spherical harmonic coefficients have m = 0. Further, 
since t{6, <fi) is even in 9, only a to where I is even are 
non-zero. This is interesting because it says that if 
the anisotropy in r is axisymmetric in the dipole frame, 
jiscr _ yrscr wnere on \y £ — multipole couplings would 
survive. While in the real world, any local free electron 
distribution is unlikely to be anisotropic in a axisym- 
metric way in the CMB dipole frame, we can see that 
in general there will be an asymmetry between odd and 
even multipoles in the power spectrum of the screening 
field. This analytic argument thus explains the "saw- 
tooth" pattern of the power spectrum of the kSZ effect 
of the Galactic disk found by Ref. [15] , 

Having developed analytic intuition for the form of 
the expected screening field, we will now present several 
physical models for the anisotropic optical depth. 
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FIG. 2: The power spectrum of the full screening field corre- 
sponding the four models for the free-electron distribution in 
the halo. The solid black curve corresponds to the general- 
ized NFW halo with a core radius of ro = 300 kpc (SphSym). 
The short-dashed blue and red curves correspond to a triax- 
ial halo (with e& — e c = 0.5) aligned with (Triax) and per- 
pendicular to (TriaxRot) the angular momentum of the disk, 
respectively. The long-dashed black curve corresponds to the 
triaxial model for the Milky Way halo proposed in Ref. [41] 
in order to explain the tidal tails of the Sgr dSph (TriaxLMJ). 
As described in Eq. |7| the amplitude of the signal depends 
linearly on the angle averaged column-density and is shown 
here for the fiducial case C = 10 21 cm~ 2 . The discussion in 
Sec. |IV| gives an analytic explanation for the sawtooth pat- 
tern found in power-spectra generated by these free-electron 
distributions. 



symmetric gaseous halo and the anisotropy is a conse- 
quence of the fact that the Solar system is offset from 
the Galactic center. The IAU standard for the Galacto- 
centric distance of the Sun is 8 kpc; however, this value 
has been recently revised upwards by ~ 5% [321 33] . We 
adopt the conservative Galactocentric distance 8.5 kpc 
to maximize the anisotropic optical depth arising from 
this offset. In the second case, it is also possible that the 
gaseous halo is triaxial, leading to an additional source 
of anisotropy. We take the virial radius of the Milky Way 
to be 300 kpc (which corresponds to an NFW halo with 
a mass of 1.5 xlO 12 M & and a concentration c = 12). 



B. Physical models for the anisotropic optical 
depth 

There are two basic ways that we may imagine an 
anisotropic optical depth through the Galaxy. In the 
first case, the Galaxy is assumed to have a spherically 



1. Spherically symmetric halo 

In order to explore the case of a spherically symmetric 
halo we need only specify the number density of electrons 
as a function of the radial distance from the Galactic 
center (GC). As discussed in Ref. [25] N-body Simula- 
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FIG. 3: Sky maps of the screening field for the three cases, 
from top to bottom: a generalized spherical NFW halo with 
a — 0.01 (SphSym); a generalized triaxial NFW halo with 
major axis perpendicular to the plane of the Galactic disk, 
a — 0.01, and = e c = 0.5 (Triax); the same as before 
but with the major axis within the plane of the Galactic disk 
(TriaxRot). In both cases the intermediate axis is chosen to 
be perpendicular to the line connecting the Solar system and 
the Galactic center. In the final case (TriaxLMJ) we take the 
model for the Galactic halo recently proposed in Ref. |41j 
in which the minor axis is aligned with the line connecting 
the Solar system and Galactic center and the major axis lies 
within the Galactic disk. 



tions which include gas dynamics generically predict a 



hot gaseous halo for Milky Way-like galaxies. For the 
spherically symmetric case we will consider a generalized 
NFW density profile with an inner core 



(1 + r/r o y 



(24) 



Spherically symmetric profiles can be written in a self- 
similar form (where the origin is at the GC) 



n e = n Q y gils (r/ro). 



(25) 



Using the Galactic coordinate system (where the origin 
coincides with the location of the Solar system and the 
GC lies on the positive x-axis) the line-of-sight density 
takes the form 

n e (t, 9, 4>) = n y gas {^Jt 2 + a 2 - 2at cos(</>) sin(0)) , 

(26) 

where t = r/r^ and a = rQc/r^. The normalization 
of the density distribution is set by requiring the mean 
angular column density be equal to 10 21 cm~ 2 . The nor- 
malization for this model can be computed analytically 

as 



1 - a 1 



n r - 



2a 2 log(a) 



2(1 -a 2 ) 2 
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2. Triaxial halos 

N-body simulations show that dark matter halos are 
not spherically symmetric but instead are better approx- 
imated by triaxial density profiles [4"4"H4"?1 . This triaxial- 
ity, if reflected in the density distribution of an extended 
gaseous halo, will produce additional anisotropy in the 
optical depth. Fits to N-body simulations [17] show that 
isodensity contours are well approximated by a radial co- 
ordinate 
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(28) 



where the cc-axis and z-axis run along the major and mi- 
nor principal axes, respectively; e& and e c (e& < e c ) are 
the ellipticities of the halo isodensity contours. Cosmo- 
logical N-body simulations give a nearly Gaussian distri- 
bution for e& with (e c ) » 0.8 and < 0.7 [UJ. We also 
note that it has been suggested the inclusion of gas cool- 
ing may make the dark-matter density distribution more 
spherical even out to the virial radius [48j . 

As with the spherically symmetric case, we consider a 
generalized triaxial NFW halo with a core. As suggested 
by Ref. [49] in order to obtain consistency between obser- 
vations of the distribution of sub-halos about the Milky 
Way and the ACDM prediction, the major axis of the 
triaxial halo must be aligned with the disk angular mo- 
mentum. In order to explore how the relationship of the 
disk angular momentum and the triaxiality of the dark 
matter halo affects the local kSZ signal, we consider the 
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two cases where the major axis is aligned with or perpen- 
dicular to the disk angular momentum. In both cases, the 
intermediate axis is along the Galactic y-axis 2 . 

We also consider the recent triaxial Milky Way halo 
proposed in Ref. |41) in order to explain the observed 
characteristics of the Sagittarius dwarf spheriodal (Sgr 
dSph). Using observations of the tidal streams of the 
Sgr dSph, Ref. [H] claim a Milky Way halo with a 
minor/major axis ratio eb ~ 0.56 and an intermedi- 
ate/major axis ratio of e c ~ 0.74, with the minor axis 
lying along the Galactic cc-axis and the major axis along 
the Galactic y-axis. 



VI. EFFECT OF HALO KSZ ON ISOTROPY 
ANOMALIES 



TABLE II: Values and p-values (computed from 10,000 
isotropic realizations) for angular momentum dispersion 
statistic Z/23 f° r the WMAP 7 year Internal Linear Combina- 
tion (ILC7) map subtracting the kinematic quadrupole (KQ) 
after accounting for the anisotropic screening from the four 
models described in the text, after scaling the maximum am- 
plitude of the screening map to 2 piK and 5 /iK. We also list 
the corresponding value for the angle-averaged free electron 
column density, C, in units of 10 21 cm -2 . 



Model 


T 1 

J-'23 


p- value 


C 




p- value 


C 




(2 /iK) 


(2 /iK) 


(2 fjK) 


(5 /iK) 


(5 /iK) 


(5 /iK) 


SphSym 


0.955 


0.0047 


12 


0.952 


0.0056 


31 


Triax 


0.962 


0.0022 


22 


0.972 


0.0008 


56 


TriaxRot 


0.950 


0.0062 


24 


0.943 


0.0091 


60 


TriaxLMJ 


0.956 


0.0042 


4 


0.958 


0.0036 
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Given a model of the anisotropic optical depth, we ob- 
tain a prediction for T scr from Eqs. ( 16 1 or ( 17 ). We can 



add that to the observed realization to obtain the actual 
cosmological realization, and calculate whatever anomaly 
statistics we want for the cosmological realization. 



A. Quadrupole— octopole planarity and alignment 

First, let us consider the angular momentum dispersion 
statistics of I 



III Taking the four models for the halo de- 



scribed in § VB labelled SphSym, Triax, TriaxRot and 



TriaxLMJ respectively, we compute the screening fields 
and add them to the WMAP ILC7 map after subtrac- 
tion of the kinematic quadrupole. In order to estimate 
the amplitude of the signal at which the p-value of the 
statistic would be affected for each case, we scale up the 
maximum of the screening map to an amplitude (1 /iK, 
2 /iK, 5 /iK) which is much smaller than the primordial 
quadrupole and octopole signals, before adding to the 
ILC7 map. The 1 \iK signal is not large enough to make 
any difference to the statistics. The results for the other 
two cases are presented for the L| 3 statistic in Table |llj 
The L\ and L\ statistics of the ILC7-KQ map are not 
anomalous before or after the kSZ correction and thus 
we do not repeat them here. 

In interpreting Table [TTJ the baseline p- value to com- 
pare with is 0.0055 for the ILC7-KQ map with no kSZ 
correction. The Table shows that for the Triax case, small 
signals at the ^2 — 5 /iK level reduce this already small 
p-value by factors of ~ 2-7! When scaled to 5 /zK, the 
TriaxRot p-value increases somewhat by factors of < 2 
and the TriaxLMJ p-value decreases by roughly the same 
factor. The SphSym screening map is unable to affect 



this statistic significantly at these signal levels. These 
tests show that a statistic such as L| 3 , which is sensitive 
to the orientation of the Galactic plane, can be highly 
sensitive to small signals oriented with this plane. 

If the maximum amplitude of the screening map is be- 
tween 2-5 fiK, the average free electron column densities 
for these maps are between 4 x 10 21 - 6 x 10 22 cm~ 2 . 
As we discuss in more detail below, though the lower 
end of this range is possible given current observational 
constraints, the upper range is unlikely 



B. Hemispherical asymmetry 

We now explore the extent to which the scattering 
of CMB photons in an extended gaseous halo can ac- 
count for the observed hemispherical asymmetry [3] . The 
statistic by which this asymmetry is inferred assumes a 
dipolar modulation of an isotropic primordial CMB tem- 
perature field, and constrains the relative amplitude of 
this modulation. The dipolar modulation leads to a cou- 
pling between modes of order i and I ± 1. As we saw in 
Sec.|IV[ the kSZ also introduces such a coupling between 
modes, so it is particularly interesting to consider how, 
through the kSZ, a screening field may impact the sig- 
nificance of the dipole modulation amplitude statistic for 
constraining hemispherical asymmetry. 

In order to do this, we will use the approach presented 
in Refs. [H2 [35]. In that work an estimator is derived 
which measures to what extent the data shows a dipolar 
modulation by looking for a non-zero coupling between 
multipoles of order i and t ± 1 up to some £ max , 



2 The Galactic coordinate system is centered on the location of 
the Solar system with the Galactic a;-axis pointing towards the 
Galactic center and the Galactic z-axis pointing in the direction 
of the angular momentum of the Galactic disk. 
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(29) 



m crc 



where 



fTT _ r-iTT , s^TT 
J l = & ' °£+l 



(30) 



a tm — 



( _ | , ; „,^; (2^ + 1) (2^ + 1) (2^ 2 + i) x 



47T 

/ V m —m 



(31) 



From the discussion in Sec. IIVI it is clear that the in- 
clusion ol Thomson scattering from free electrons in an 
extended Galactic halo will couple multipoles of order £ 
to those of order I — 1 and i + 1 in the observed ae m s. 
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FIG. 4: How a screening field can affect the measurement of a 
hemispherical asymmetry. The grey band indicates the con- 
fidence level of the distribution at that value of f ma x ranging 
from 1 a to 3 a. The solid thick black curve shows the value 
of i&i evaluated on the ILC cleaned WMAP7 map. From top 
to bottom we show the resulting hemispherical asymmetry for 
TriaxLMJ, SphSym, TriaxRot, and Triax as a function of £ ma x. 
This Figure demonstrates how the orientation of the triaxial- 
ity of the halo greatly affects its ability to impact the inferred 
amplitude of a hemispherical asymmetry. 

To check the effect of the kSZ on the significance of 
the hemispherical asymmetry, we evaluate the estimator 
wf T on the inferred cosmological signal given by Tf™ = 
Tim + ^Tm(^) f° r various values of the angle-averaged 
electron column density C. As in the previous section, the 
four extended halos considered here are the spherically 
symmetric case (SphSym), and the three triaxial halos 
with various orientations (Triax, TriaxRot, TriaxLMJ). 

We show the value of the estimator w\ evaluated on 
the WMAP ILC7 map in Fig. B] Because of the scale 
dependence of the power spectra of the screening fields, 
seen in Fig. [51 we plot the value of wi as a function of 
f max in Fig. Hi The grey band indicates the confidence 



level of the distribution at that value of £ max ranging 
from 1 a to 3 cr. The Figure shows the value of w\ for 
(from top to bottom) for TriaxLMJ, SphSym, TriaxRot, 
and Triax. It is clear that the orientation of any triaxi- 
ality has a significant impact on the inferred amplitude 
of the hemispherical asymmetry. In particular, the Tri- 
axLMJ orientation leads to an increase in the inferred 
amplitude. This can be seen in Fig. [3] since the inclu- 
sion of the TriaxLMJ screening field adds power in the 
northern Galactic hemisphere and subtracts it from the 
southern hemisphere. For the other cases, the inferred 
amplitude decreases. 

Since the anisotropies induced by a local kSZ are at the 
level of 1CT 5 fj,K for I > 10 (see Fig. [2|, in order for the 
screening field to have a significant impact on the inferred 
amplitude of w\, we must have C/(10 21 cm~ 2 ) ~ 10 4 . A 
free-electron column density of this magnitude is ruled 
out by observations of OVII and OVIII absorption [28] 
and pulsar observations towards the Large Magellanic 
Cloud (LMC) !29j. Therefore, we find that a local kSZ 
is unlikely to provide a plausible explanation for a hemi- 
spherical asymmetry. 



VII. DISCUSSION 

Although measurements of the CMB have generally 
confirmed our current understanding of the formation 
and evolution of the universe, there are several anoma- 
lies which appear to be at odds with our standard pic- 
ture. These anomalies could suggest that there may be 
some process which violates the statistical isotropy of the 
CMB fluctuations on large scales. Although attempts 
have been made to modify the physics of the early uni- 
verse in order to explain these anomalies, there are also 
strong reasons to look for explanations in the local uni- 
verse. First, the directionality of these anomalies seems 
to closely coincide with both the Solar CMB dipole as 
well as the velocity of the Galactic barycenter relative 
to the CMB rest-frame. Second, given that the anoma- 
lies are on large angular scales, any non-primordial and 
causal explanation must be local in origin. 

Here we have explored how a local kSZ signal (due to 
the motion of any extended hot gaseous halo associated 
with the Milky Way relative to the CMB) may affect the 
significance of several anomalies observed in the CMB 
anisotropies. Both theoretical and observational con- 
siderations indicate the presence of a hot gaseous halo 
with an extent of several tens of kiloparsecs [551 ED]- 
Anisotropies in the optical depth through this gaseous 
halo can be due to the offset of the Solar system from the 
Galactic center, as well as any triaxiality in its distribu- 
tion. We computed the kSZ signal from several plausible 
physical models for the shape and orientation of the halo, 
and studied their impact on the observed CMB sky, in 
the form of the WMAP ILC7 map. 

We considered how a local kSZ may affect the observed 
planarity of the CMB quadrupole and octopole moments 
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as well as their relative alignment, motivated by the fact 
that kSZ from a triaxial halo can naturally relate to the 
special directions that appear to be associated with the 
£ = 2 and 3 multipole moments. Surprisingly, we found 
that relatively small changes in the observed amplitude 
of these moments (~ 10%) to account for a kSZ signal 
can reduce the already tiny p-values for these anomalies 
of up to factors of 2-7. The corresponding free electron 
column density needed to affect the planarity/alignment 
statistics at this level is 4 - 60 x 10 21 cm~ 2 . 

A local kSZ signal couples multipole moments of or- 
der I to those of order £ ± 1, so it is natural to consider 
how such a signal would affect any inferred hemispher- 
ical asymmetry in the CMB anisotropics in the form of 
a dipolar modulation, which gives rise to the same cou- 
pling. Using a statistic first derived in Ref. [TS] , we found 
that a local kSZ signal can have a significant effect on the 
inferred amplitude of a dipolar modulation of the primary 
CMB anisotropies. In this case, the corresponding free 
electron column density would need to be ~ 10 25 cm~ 2 . 

Theoretical and observational constraints on the free 
electron fraction in a extended hot gaseous halo asso- 
ciated with the Milky Way place a fairly strict bound 
on the column density of < 10 21 cm" 2 . The most pre- 
cise constraints come from observations of the dispersion 
measure to individual pulsars |29) which find a column 
density in free electrons to the LMC of ~ 3 x 10 20 cm~ 2 . 
Given that the LMC is ~ 50 kpc from the Galactic center 
and that the scale-radius of the Milky Way halo is ~ 300 
kpc, these observations indicate that the total optical 
depth may be as large as 10 21 cm -2 . In addition to this, 
observations of local (z = 0) OVII and OVIII absorption 
towards several quasars indicates the presence of an ex- 
tended hot gaseous halo around the Milky Way [10] . The 
inferred free electron column density is quite uncertain 
given assumptions about the metallicity of the gas (Solar 
abundance was assumed to convert the observations to 
an electron density; sub-Solar values would increase the 
inferred density), as well as a model dependence coming 
from the assumed profile of the gas. Finally, theoreti- 
cal considerations imply the existence of a hot extended 



gaseous halo with a fractional mass of the order of the 
cosmic baryon fraction, / ~ 0.1. Therefore, for a Milky 
Way sized halo (M ~ 1.5 x 10 12 M Q ) with a scale ra- 
dius of 300 kpc, we would expect a column density no 
larger than 10 21 cm -2 . Given that, for a local kSZ signal 
to have a significant effect on the CMB anisotropies, we 
must have a free electron column density > 10 21 cm' 2 , 
but that both theoretical and observational considera- 
tions place the limit at < 10 21 cm -2 , it is unlikely that a 
local kSZ signal can explain any of the CMB anomalies 
considered here. Even given the uncertainties, the kSZ 
signal can at best only be at the lowest amplitude needed 
to affect the anomaly statistics. 

Finally, we note that contamination from kSZ in the 
Solar system can also, in principle, provide an anisotropic 
contamination of the CMB sky at large angles. The ge- 
ometry of the heliopause is coincidentally aligned with 
the CMB dipole [TB], providing a motivation for looking 
there for a explanation for the special directional prop- 
erties of the large-angle isotropy anomalies. However, in 
practice, the optical depth of free electrons in the Solar 
system is more than 7 orders of magnitude below what 
is required to produce the necessary signal. 
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